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ABSTRACT 

Gas falling quasi- spherically onto a black hole forms an inner accretion disc if its specific 
angular momentum Z exceeds Z» ~ r g c where r g is the Schwarzschild radius. The standard 
disc model assumes I 3> l*. We argue that, in many black-hole sources, accretion flows have 
angular momenta just above the threshold for disc formation, I <; Z„, and assess the accretion 
mechanism in this regime. In a range I* < I < Z cr , a small-scale disc forms in which gas spirals 
fast into the black hole without any help of horizontal viscous stresses. Such an 'inviscid' disc, 
however, interacts inelastically with the feeding infall. The disc-infall interaction determines 
the dynamics and luminosity of the accretion flow. The inviscid disc radius can be as large as 
14r s , and the energy release peaks at 2r g . The disc emits a Comptonized X-ray spectrum with 
a break at ~ 100 keV. This accretion regime is likely to take place in wind-fed X-ray binaries 
and is also possible in active galactic nuclei. 

Key words: accretion, accretion discs — binaries: general — black hole physics — radiative 
mechanisms: thermal — X-rays: galaxies — X-rays: stars 



1 INTRODUCTION 

Spherical accretion flows onto black holes (BHs) have low ra- 
diative efficiencies (Shvartsman 1971; Shapiro 1973; but see also 
Bisnovatyi-Kogan & Ruzmaikin 1974; Meszaros 1975). If the flow 
is slightly rotating, the situation changes crucially: a small-scale 
accretion disc forms and then the BH 'switches on' as a luminous 
source. The condition for disc formation reads Z > Z* where Z is 
the specific angular momentum of the accreting matter and Z* = 
0.75r 9 c (see Section 2.2), r g — 2GM/c? being the Schwarzschild 
radius. The standard disc model assumes I ^> Z» (see Pringle 1981 
for a review). High-Z discs can form, e.g., in X-ray binaries fed by 
streams from Lagrangian L\ point. 

By contrast, in wind-fed X-ray binaries, the accretion flows 
are quasi-spherical with Z ~ Z» (see e.g. Frank, King, & Raine 
1992 for a review). The average Z of the accreting wind matter is 
1 = f(l/4)£lR^ (Illarionov & Sunyaev 1975; Shapiro & Lightman 
1976) where R a = 2GM/w 2 ~ 10 11 cm is the accretion radius, 
w w 10 8 cm s" 1 is the wind velocity, Q is the angular velocity of 
the binary, and ( ~ 1 is a numerical factor. Three X-ray binaries 
have been classified as massive ones with black hole companions: 
Cyg X-l, LMC X-l, and LMC X-3 (see Tanaka & Lewin 1995). 
They have orbital periods P = 2n/Q, — 5.6, 4.2, and 1.7 d, re- 
spectively. Hence, T 1.5CP" 1 (A//M )O/lO 8 )" 4 r 9 c ~ L in 
these objects, that implies marginal formation of a small-scale disc. 
Whether the disc forms or not depends on the precise value of an- 
gular momentum which is difficult to calculate because the details 
of wind accretion are not fully understood (note the strong depen- 
dence of Z on the wind velocity, w). The fact that the objects are 
luminous indicates that a disc does form. 



Small-scale accretion discs may also form around massive 
black holes in active galactic nuclei (AGN). Among possible 
sources of the accreting gas in AGN are star-star collisions and 
tidal disruption of stars by the BH (Hills 1975). The accretion flow 
is then likely to have a modest angular momentum. 

The regime Z <; Z* can be widespread among luminous black- 
hole sources. Note that only a small number of wind-fed BH bina- 
ries and galactic nuclei are active, and most of the unseen objects 
may accrete in the spherical regime Z < Z» with a low luminosity. 
The observed bright sources with Z > Z» are then in the 'tail' of 
the distribution of objects over angular momentum, dN/dl. If this 
distribution falls steeply towards high Z, most of the bright sources 
should be near the threshold for disc formation. 

In this paper, we construct an axisymmetric model of accretion 
with Z fZ, Z„. The outline of the model is as follows. The accretion 
flow starts at large distances where it is almost spherical; the gas is 
Compton cooled and it falls freely towards the BH. Near the BH, 
the rotation deflects the infall from the radial direction and a disc- 
like caustic appears in the defocused flow. Here the flow liberates 
energy in inelastic collision and then proceeds via the thin disc into 
the BH. Such a disc is drastically different from its standard coun- 
terpart as regards the dynamics, the mechanism of energy release 
and the emission spectrum. 

The paper is organised as follows. In Section 2, we study the 
axisymmetric free fall in the gravitational field of a Schwarzschild 
BH and discuss the caustic in the equatorial plane (the collision 
ring) where the accretion disc is formed. In Section 3, we write 
down the basic equations of the disc. In Section 4, we solve the 
equations numerically, obtain a model for the disc structure, and 
find the efficiency of accretion. Most of the energy is released quite 
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close to the event horizon, at ~ 2r g , and the problem requires exact 
relativistic treatment. All the calculations are relativistic (assuming 
a non-rotating BH). For comparison, in Section 5, we also formu- 
late and compute the analogous Newtonian problem. The results 
and issues for further study are discussed in Section 6. 



2 FORMATION OF A SMALL-SCALE DISC 

Let us consider an axisymmetric accretion flow slowly rotating 
around the polar axis, so that the non-radial velocity v v <S v T at 
R 3> r g (an asymptotically spherical inflow). As a result of Comp- 
ton cooling by the inner source, the inflow becomes super-sonic 
inside the Compton radius, Rc ~ 10 4 r 9 (see e.g. Illarionov & 
Kompaneets 1990). AtR<Rc, the matter accretes along ballistic 
parabolic trajectories. The freely falling flow in Newtonian grav- 
ity is discussed in Illarionov & Beloborodov (2000, hereafter IB). 
Here we need the relativistic equations since the typical radius at 
which accretion switches from the spherical to a disc-like regime, 
R ~ I 2 /r g c 2 , is comparable to r g in our problem. Hereafter label 
'oo' corresponds to radii R^> r g (but R < Rc). 



2.1 Free fall 

Consider a free-fall trajectory in Schwarzschild (t, R, 9, ip)- 
coordinates. Let A be proper time along the trajectory. Denote the 
orbital angular momentum by 1; its projection onto the polar axis is 

l z = / sin #00. 

A parabolic trajectory (with specific orbital energy E — c 2 ) is de- 
scribed by equations (see e.g. Misner, Thorne & Weeler 1973) 



J/3 



ctg 2 0, 



R 



• dip 
d\ 



(1) 



(2) 



(3) 



The only difference from the corresponding Newtonian equations 
is the presence of the additional factor (1 — r g /R) in equation (1). 
One can solve equations (1-3) for the trajectory 9(R), <p(R), 



cos 8 — cos #oo cos ip, 
sin tp 



sin(ip - ipoo) 



HR) = 



sin 9 ' 



IdR 



sJr g c 2 R 3 -l 2 R 2 (l-r g /R) 



(4) 
(5) 

(6) 



The ip has the meaning of an orbital angle. In the Newtonian limit, 
r g /R — > 0, equation (6) gives cosi/> = 1 — l 2 /GMR (eq. 11 in 
IB). 

The accreting matter streams from a sphere Soo along the bal- 
listic trajectories. The flow is assumed to be symmetric about the 
equatorial plane. In this plane, a streamline has to collide with its 
symmetric counterpart (see Fig. 1). The loci of collisions form a 
two-dimensional caustic, the collision ring. 

We assume that there are no intersections of the sreamlines 
before they reach the equatorial plane. It is so if the Jacobian of 
the transformation flooj^oo — > 0(R),ip(R) stays positive along 
the streamlines. Taking into account that dOii/dipoo = and 




Figure 1. Schematic picture of the disc formation. The inflow through 
Soo has angular momentum I increasing towards the equatorial plane, 
dl/d sin #00 ^> 0. The collision radius, which is a growing function of 
I, monotonically increases as Poo approaches n/2. The shadowed parts of 
the accretion flow collide outside r g and form a couple of radiative shocks 
which sandwich the collision ring. The other part of the flow (at small polar 
angles where I < U) is directly swallowed by the black hole with a low 
radiative efficiency. 



d<pn/d<poo — 1 for any axisymmetric flow, the condition of posi- 
tive Jacobian reads 



d cos 9 R /j ■ / dtp 

= cos ip — cos 9oo sin ip — > 0. 



d cos 9 



d cos 9 



It is satisfied if cos# 00 dV'/dcos6 l 00 < 0. Since ip increases with 
increasing I, the condition can be rewritten as 



dl 



> 0. 



(7) 



d sin 9 a 

Hereafter we assume that condition (7) is satisfied. 
2.2 The collision ring 

Streamlines coming from Soo with 9^ k/2 collide in the sym- 
metry plane and define a mapping foo , 9^ — > r, ip. Here (r, ip) 
are polar coordinates in the equatorial plane. (The r has the same 
meaning as R; we use different notation to emphasize that r mea- 
sures radius in the plane.) The mapping is given by equations (4-6) 
with 9 = 7r/2. Equations (4) and (5) say that a streamline reaches 
the equatorial plane at ip = 7r/2 and executes 1/4 of the full turn 
around the polar axis, ip = (poc + n/2. A circle 9oo = const at 
Soo transforms into a circle r = const with the 7r/2 rotation. The 
radius of collision is found from ip = ty/2, 



IdR 



^Jr g c 2 R 3 -PR 2 (l-r g /R) 



(8) 



Here / (#00 ) = l(n — 9^ ) is the angular momentum of the colliding 
streamlines. Equation (8) involves an elliptic integral. At r r g , it 
is simplified and yields the Newtonian relation, r — I / GM. The 
relativistic effects do not crucially alter the Newtonian expression 
even at r ~ r g . The solution to equation (8) can be written as an 
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expansion which, in the leading order, yields 



I 2 (10-3tt) 



I 2 
GM 



GM 4 " GM 7 

This equation reproduces the exact solution r(l) with high accuracy 
(the maximum error ^ 0.5 per cent at r — r g ). 

The collision occurs outside the event horizon if I exceeds 
a minimum value, h, such that r(U) = r g . From equation (8) 
one finds I, ^ 0.754r 9 c. (In the Newtonian approximation, h = 
r g c/V2 « 0.707r 9 c.) A streamline with I < I, goes directly 
into the BH (see Fig. 1). A streamline with 0.754 < l/r g c < 2 
would also end up in the BH (see e.g. Novikov & Frolov 1989). 
The collision in the equatorial plane, however, happens before that, 
at r <, 8r g . 

Streamlines starting at 9cc — » 7r/2 have the largest angular 
momentum, I — » lo, and define the outer edge of the collision ring, 
r « ll/GM - r g /7. 

A collisionless shock forms on each side of the ring. The ac- 
cretion picture shown in Fig. 1 assumes that the shock is radiative 
and pinned to the ring; we study this condition in Section 2.3. The 
free-fall equations (1-3) give the three components of the flow ve- 
locity in the preshock region; note that the azimuthal velocity is 
about the radial one. 



2.3 The pinned shocks and vertical mixing 

The vertical component of the plasma velocity is reduced at the 
shock front by a factor of £ > 4, while the horizontal velocity re- 
mains unchanged (the horizontal motion towards the BH can stay 
super-sonic). The dissipated kinetic energy goes mostly to ion heat- 
ing and some fraction e may go to electron heating (studies of inter- 
planetary and interstellar shocks indicate e ~ 0.2 — 0.3, see Draine 
& McKee 1993). The electrons immediately (at the shock front) ra- 
diate the received energy via inverse Compton scattering, produc- 
ing hard gamma-rays (Shapiro & Salpeter 1975). The gamma-rays 
are converted into e ± pairs and increase the downstream electron 
density by a factor of z <; 1. The cooled downstream e ± then take 
the energy from the ions. 

The radiative capability of the shocked plasma depends on e 
and also on the rate of energy exchange between the downstream 
electrons and protons. The minimal e-p coupling is provided by 
Coulomb collisions: the electrons of temperature T e cool the pro- 
tons of temperature T p on a time-scale (Landau & Lifshitz 1981), 

- I n/2m p 9^ 2 



tCoul 



.T, 



3/2 



17^ 



-s, 



(9) 



m e ccjTzn In A zn 

where n is the proton density, and a Coulomb logarithm In A = 15 
and e = kT e /m e c 2 < 1 are assumed. The e-p coupling can be 
enhanced by collective plasma effects, so that the proton cooling 
time-scale is t cp = x _1 ^Coui where \ > 1. 

The shock is radiative and held-down to the equatorial plane 
if t cp is shorter than £,r/vs where vs = c(r 9 /r) 1 ' 2 is the typ- 
ical free-fall velocity. The downstream density can be estimated 
in terms of the accretion rate M to t- n w £M to t /^r 2 m p vs. The 
condition that the shock is pinned to the equatorial plane then reads 



2imi p 8\ 



3/2 



zln Am e £ 2 xrh 



< 1, 



(10) 



where rh = M to tc 2 /Lb is the dimensionless total accretion rate 



and Lb = 27rr £ 



/<tt is the Eddington luminosity. 



The downstream electron temperature can be found from the 
energy balance, 



zTe 

tc 



te 



(ID 



where tc = 3m e c/8a T U is the Compton cooling time-scale 
and U is the radiation energy density. For a radiative shock, 
the proton thermal energy is transformed into a radiation flux 
F ~ (3/2)nkT p v s /£. It yields an estimate, U ~ 2F/c ~ 
3nkT p VB /c£, and hence 



tc 



m e c 2 £, 



8a T nkT p Vf{ ' 
The energy balance (11) then yields 

2/5 



lnAm e £x c 
2iTm p 4«ff 



0,1 (4w ff ) 



,\ 2/5 



(12) 



we oc r 



-1/2 



and hence e oc r 1 ^ 5 . The temperature weakly de- 
pends on parameters. 

The hot downstream layer radiating most of the proton heat 
has optical depth tt ~ (vs / ' £)t e -pZnaT . It satisfies the relation, 
rrOe ~ 0.1 (as follows from eqs. 9 and 12), i.e, tt ~ 1 for typ- 
ical 9 e ~ 0.1. With more exact approach, the unsaturated Comp- 
tonization in the hot layer involves radiative transfer with multiple 
scatterings, and tt ~ 1 corresponds to a Kompaneets' parameter 
y ~ 1 (see e.g. Rybicki & Lightman 1979). In deeper layers, at 
tt > 1, the downstream temperature (both T p and T e ) falls off 
sharply (see Zel'dovich & Shakura 1969, Shapiro & Salpeter 1975 
for a similar shock on the surface of a neutron star). 

Using equation (12), we rewrite the condition (10) as 



rh > rho 



0.5 

zx 2/s U 



7/5 



(±Y" > „™*(i) (±) .as) 

V va J z9 e \ £ J \vttJ 



We are interested here in sub-Eddington sources (with luminosities 
L < Le) so that the free-fall approximation can be used for the 
super-sonic infall. It requires rh < ibe = n^ 1 where n is the 
radiative efficiency of the disc (see Section 4 and Fig. 6). In the 
range rho < rh < ij" 1 we have sub-Eddington accretion with 
radiatively efficient shock pinned to the equatorial plane. 

This pattern of accretion is different from the previously con- 
sidered situation rh ^ rho (see e.g. Igumenshchev, Illarionov & 
Abramowicz 1999). In that case, the plasma is unable to cool fast 
and supports quasi-spherical shocks away from the collision ring. 
The hydrodynamic simulations of Igumenshchev et al. (1999) show 
that the shock tends to shrink towards the ring of radius l 2 /GM 
when rh increases. 

We hereafter assume rh > rho and have a free-falling flow 
until it reaches the thin disc of shocked plasma. The disc has a finite 
half-thickness H supported by the pressure pd ~ PdC 2 against the 
ram pressure of the infall, p ft ~ PsVg. Here p d and c s are the 
density and the sound speed in the midplane of the disc, and ps 
is the density of the infall. The pressure balance yields pd/ ps ~ 
(vf{/c s ) 2 . The regime of pinned shock implies efficient radiative 
cooling and we have c s vg and p d ^S> pa- The accretion velocity 
in the disc is comparable to vs (see Sections 3 and 4) and from 
matter conservation we estimate H/r ~ ps/ 1 p d ~ (c s /vs ) 2 . 

At a given radius, the disc is composed of matter that en- 
tered the disc at larger radii with different horizontal velocities. A 
strong turbulence is likely to develop in the disc under such condi- 
tions. The turbulence can mix the disc in the vertical direction on 
a time-scale ~ H/c s which is ~ c s /vs shorter than the accretion 
time-scale, t a ~ r/vg. In the mixing process, the random motions 
above the level ~ pdC 3 are converted into heat and the disc acquires 
a vertically-averaged horizontal velocity (with both radial and az- 
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imuthal components). The released heat is radiated away, contribut- 
ing to the total luminosity. 

When absorbing the Mailing matter, the disc also absorbs hor- 
izontal momentum and its velocity changes. The process of disc- 
infall inelastic collision is governed by the corresponding law of 
momentum conservation (see Section 3). The disc is a 'sticky' caus- 
tic in the accretion flow and it resembles the accretion line in the 
Bondi-Hoyle-Lyttleton (BHL) problem (see Bondi & Hoyle 1944). 
We discuss this analogy in more detail below. 



3 FAST INVISCID ACCRETION DISC 

The shocked turbulent gas continues to accrete into the BH through 
the thin disc. We will show that if the inflow angular momentum is 
smaller than some specific value then accretion in the disc occurs 
fast: the gas spirals into the BH on the free-fall time-scale. In con- 
trast to the standard viscous disc, the radial velocity is much larger 
than the sound speed. Therefore, turbulent viscosity has a negli- 
gible effect on accretion, i.e. the disc is inviscid in the horizontal 
directions (yet the turbulence efficiently mixes the disc in the verti- 
cal direction, see Section 2.3). Thereafter we solve the problem of 
steady disc accretion in the inviscid regime. In illustrative calcula- 
tions we assume solid body rotation of the inflow at Soo , 

t/2 2 

l{0oo) = lo sinfloo, h = I I sin 2 OoodOoo = (14) 



/>7T 

Jo 



Here l z is the average specific angular momentum of the inflow. 



3.1 Accretion rate 

The asymptotically spherical inflow at Soo has a total accretion 
rate M to t and a homogeneous distribution dM/d^loo = M to t/47r 
where dO,oo = d^ood cos 9oo . All the matter impinging the disc fi- 
nally accretes into the BH. Hence, the disc accretion rate at a radius 
r, M(r), equals the rate of matter supply to the disc at all r' > r, 
which in turn equals the rate of accretion with I' > / through the 
equatorial segment sin 8'^ > sin 8oo at Soo (see Fig. 1), 



M(r) = M >t = M tot cos9oo{l). 



(15) 



In the case of solid body rotation at Soc (eq. 14), we have the 
relation sin 9oo = l/lo and hence 



(16) 



Here I is a function of r, see equation (8). With the approximate 
relation I 2 « GMr, we get 



M(r) «Mtot «/l 



ro 



(17) 



3.2 Angular momentum 

The specific angular momentum of a steady inviscid disc, j(r), 
equals the average angular momentum, l' z , of matter absorbed by 
the disc at r' > r. (It follows from the conservation of the accret- 
ing mass and conservation of the total accreted angular momentum. 
The radiation carries away a small fraction of the disc momentum, 
see Appendix A; hereafter we neglect these losses.) We get, 

M{r) J r dr' M>; ] l dV 



Here l' z — V sin 6^ and M>; is given by equation (15). The infall 
has lz(r) < j(r) and decelerates the disc rotation: j(r) decreases 
inwards. 

In the case of solid body rotation at Soo, we have sin — 
l/l , h = l 2 /lo and dM/dl = -M tot (l/l )(ll ~ l 2 )~ 1/2 - Work- 
ing out the integral (18) then yields 

9 1 9 l 2 (r\ 

iW = ^o + i/ z M = ^o + ^, (19) 

(see Fig. 2). With the approximate relation / w GMr, we get 



(20) 



3.3 Radial equation 

The matter in the disc is not in free fall because it absorbs the in- 
falling matter. The infall accelerates the accretion in the disc (the 
law of momentum conservation is given in Appendix A). The equa- 
tion for the disc radial velocity u r reads, 



du r 
dr 



dM (u r - u r ) r 3 c 2 



dr 



M{r) 



2r 2 u r 



3l 



Here j is given by equation (18) and 

2 2 

2 rr g c 
JK 2r - 3r 



(21) 



(22) 



At r > (3/2)r 3 , j'k nas the meaning of the squared angular mo- 
mentum for circular Keplerian rotation. At r < (3/2)r s , such mo- 
tion is not possible in Schwarzschild geometry, which is reflected 
by the fact that j K < 0. The u r in equation (21) is the radial veloc- 
ity, dr/dX, of matter impinging the disc at radius r (see eq. 1), 



1 • 



(23) 



with I being related to r via equation (8). 

Note that equation (21) is formally identical to the correspond- 
ing equation in Newtonian gravity (considering u r as a normal 
Newtonian radial velocity). The relativistic effects come into the 
problem through j\, u r , and the relation (8) between I and r. 

The first term on the right-hand-side of the radial equation is 
a result of the disc-infall interaction. Assume for a moment that 
inside a radius n there is no matter supply to the disc. Then at 
r < n the disc matter is freely falling with constant M, j, and E, 
where E is the orbital energy given by equation (24) below. In the 
presence of matter supply, dM /dr 7^ 0, both j and E change with 
radius. 

The radial equation (21) is similar to the BHL equation de- 
scribing gas motion along the accretion line (see eq. 8 in Bondi 
& Hoyle 1944). The only difference is that we have an additional 
term proportional to j 2 j This term reflects the fact that the disc 
matter is rotating and a repulsive centrifugal force appears. The 
centrifugal force dominates over gravity if j > j'k and it can stop 
accretion. Therefore, a steady inviscid disc may form only in accre- 
tion flows with limited net angular momentum (see Section 3.5). 

The differential equation (21), like the BHL equation, does 
not admit analytical solution. We solve numerically equation (21) 
combined with the formulae (15), (18), (22), (23) for M, j, 
u r ', respectively, and account for the relation (8) between I and r. 
The outer boundary condition is u r (ro) = ti~ . The solution u r (r), 
j{r), M(r) yield a closed description of the fast inviscid disc. The 
model is fully determined by the inflow angular momentum l{6oo)- 
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3.4 Energy release 



Once the disc dynamics is solved and u r (r), j(r), M(r) are 
known, one can compute the energy liberated by the disc-infall col- 
lision. The orbital energy of the disc matter E — —u t (see Appen- 
dix A) equals c 2 at the outer edge of the disc and decreases inwards. 



From UiU 1 



-c we have 



E 2 = c 2 (u r f + V 2 , 



V 2 (r) 



<? + h 



(24) 



The material absorbed by the fast disc at r' > r and accreted down 
to r liberates a luminosity L >r — M(r)[c 2 — E(r)]. The total 
energy released by the disc is 

2 



L = M d [c 2 - E{r g )] 



(25) 



Here Md = M(r g ) is the total accretion rate through the disc and 
c 2 — E(r g ) is the binding energy of gas swallowed by the BH. The 
energy liberated per unit radius is 



dr 



-[M(r) (c 2 -E(r))] 



(26) 



3.5 The limitation on the inflow angular momentum 

The essential condition for an inviscid regime of accretion is that 
there is no turning point in the disc (u r < everywhere). It re- 
quires E > V (see eq. 24) where V(r) is the effective potential 
for radial motion (Landau & Lifshitz 1975; see Misner et al. 1973, 
ch. 25 for detailed discussion). The potential has a pit at small r and 
matter may fall into the BH. With increasing angular momentum, 
there appears a centrifugal barrier enclosing the pit and preventing 
accretion. 

The main quantity determining the disc dynamics is the angu- 
lar momentum j(r). The no-turning-point condition E > V can be 
expressed as 



j(r) < j t (r) ee 




(27) 



where j = j t corresponds to E = V . The disc orbital energy is not 
much below the parabolic value c 2 and, in the first approximation, 
jt can be estimated assuming E = c 2 , 



jto = crJ—2—, (28) 



at r = 2r g (see Fig. 2). The fast disc can exist only with j(2r g ) < 
2r g c. Note that j(r) > l z at any r, where l z is the average angular 
momentum of the inflow. Therefore, the necessary condition for the 
fast-disc regime reads 



l z < 2r g c. 



(29) 



This limitation does not exclude distributions 1(8 oc) with a tail of 
large I at 6 X — > 7r/2 (and the inviscid disc can have an infinitely 
large radius). However, most of the accreting mass is required to 
have a modest angular momentum. 

If the condition (29) is not satisfied, there is no stationary so- 
lution to the inviscid problem. Instead a viscous disc forms which 
spreads out to large radii (Kolykhalov & Sunyaev 1979). The tran- 
sition to the viscous regime is rather complicated, showing unstable 
time-dependent behaviour, and it will be studied in a separate pa- 
per (Beloborodov & Illarionov, in preparation). Here we restrict our 
consideration to the fast inviscid regime. 




Figure 2. Two characteristic angular momenta of the problem, jto(r) 
(eq. 28) and juir) (eq. 22), are shown by the heavy and dashed curves, 
respectively. The necessary condition for the fast disc regime is that the 
disc angular momentum j(r) < Jto(r) everywhere. The bottleneck for 
accretion is at 2r g . The solid lines show j(r) in the model computed in 
Section 4 for two cases: io = 'cr ~ 2.62r 9 c (marked 1) and in = 2r 9 c 
(marked 2). The angular momentum of matter impinging the disc, l z (r), is 
shown by dotted lines. 



4 NUMERICAL MODEL 

We now compute the illustrative model of a fast disc formed by the 
inflow with angular momentum (14); then Iq = (3/2)T z is the only 
parameter of the problem. The limitation (29) requires lo < 3r 9 c. 
The exact critical value of l is found below. 



4.1 Dynamics 

The disc accretion rate M(r) and its angular momentum j(r) are 
given by equations (16) and (19). 

The solution for the accretion velocity is shown in Fig. 3. At 
l — l CI m 2.62r 9 c, the velocity profile touches zero: accretion is 
stopped by the centrifugal barrier at r cr w 2.14r 9 (j touches j t , 
see eq. 27, and it always stays smaller than jto, see Fig. 2). The 
critical radius is near 2r 9 , close to the minimum of jto- At I > l cr 
the fast inviscid disc cannot exist and hence its maximum radius is 
r max = ro(Zcr) w 13.6r 9 « 27GM/c 2 . 

In Fig. 4, we show the trajectory of the spiraling gas in the 
disc. In the critical case, the gas makes infinite number of revo- 
lutions at r C i before it falls into the BH. The disc surface density 
E(r) = M(r)/2nru r — > oo at r — > r cr . 



4.2 The radiative efficiency 

In Fig. 5 we show the distribution of the liberated energy over ra- 
dius, rdL/dr. The energy release peaks at w 2r 9 . The sharp peak 
in the critical case is caused by the minimum of u r (see Fig. 3). The 
relative velocity between the disc and the infall increases dramat- 
ically at r = r cr and the inelastic collision liberates more energy. 
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Figure 3. The radial velocity of the disc, u r (in units of c). The long- 
dashed, dashed, and dotted curves correspond to the cases lo = la ~ 
2.62r 9 c, (o = 2r 9 c, and (rj = \f2r g c, respectively (see eq. 14 for the 
definition of lo). The outer radius of the disc, rrj ~ 2l^/r g c 2 — f g /7, is 
shown by solid dots. The radial velocity of matter impinging the disc, u r , 
is shown by the solid curve (eq. 23). 



Figure 5. The radial distribution of the disc luminosity: solid curves - 
the released luminosity (eq. 26), dotted curves - the observed luminosity 
corrected for the gravitational capture by the black hole (eq. 30). Labels 1, 
2, and 3 correspond to the cases irj = l CI ra 2.62r 9 c, lo = 2r 9 c, and 
lo = \/2r g c, respectively. 



-i — i — I — i — i — r~ 



1 =2'/V 



I I I I I I I I I I i l l 




Note that in the standard disc model (Novikov & Thorne 1973) 
rdL/dr peaks at ~ 9r 9 . 

The small size of the main emitting region implies that the 
released energy is partly captured into the BH. Accounting for the 
capture effect, we evaluate the observed luminosity in Appendix B, 



di 



dr 



n(r)} . 



(30) 



The capture is enhanced by the Doppler beaming of the disc radi- 
ation into the BH: the emitting material accretes fast. The effect is 
especially strong when lo is small (see Fig. 5). In the critical case, 
lo = Icr, the accretion is decelerated at r w 2r 9 and the luminosity 
capture is less efficient. 

There are a few ways to define the radiative efficiency of 
accretion: (i) r\ = L/M tot c 2 , (ii) r\ = L/MdC 2 , (iii) 77 = 
Lobs/MtotC 2 , and (iv) 77 = L a b s /MdC 2 . We choose the first def- 
inition as a basic one and compute the other three for comparison. 
The results are shown in Fig. 6. 

The maximum efficiency, r\ w 0.0372, is achieved in the crit- 
ical case, lo — l C r ■ It is smaller than the efficiency of the standard 
disc around a Schwarzschild black hole, r\ = 1 — 2\/2/3 w 0.0572, 
however it is comparable to this value. We conclude that the fast in- 
viscid disc is quite an efficient regime of accretion. 



Figure 4. The trajectory of gas motion in the disc, from the outer edge into 
the black hole. In the critical case (bottom panel) the trajectory makes an 
infinite number of revolutions at the critical radius r cr pb 2.14r 9 before 
falling into the black hole. 



5 NEWTONIAN FAST DISC 

In this section, we formulate and solve the fast disc problem in 
Newtonian gravity; the importance of the relativistic effects is well 
seen when comparing the BH fast disc with its Newtonian coun- 
terpart. The natural units of length and velocity in the relativistic 
problem are r g and c. The corresponding units in the Newtonian 
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Figure 6. The efficiency of accretion, r\, versus l . The four curves corre- 
spond to different definitions of efficiency: solid curve - r) = L/Mtotc 2 , 
dashed curve - r) = L/M d c 2 = 1 — E(r g )/c 2 , dotted curve - r) = 
L Qas / Mtotc 2 , l° n g dashed curve — rj = L oos / M d c 2 . Here L is the power 
released by the disc (eq. 25), L bs is the luminosity corrected for the cap- 
ture by the black hole, JWtot is the total accretion rate (including the polar 
inflow that plunges into the event horizon before crossing the equatorial 
plane, see Fig. 1), and M d is the accretion rate through the disc. 



problem are the accretor radius R r and the free-fall velocity 

/2GM\ 1/2 



V R* ) 



5.1 Dynamics and energy release 

The Newtonian disc fed by matter coming from Soc with angular 
momentum l(Qoo) is described by the same dynamical equations 
(15), (18), and (21) as the relativistic disc. The only difference ap- 
pears in the expressions for r(l), j^(r), and u T , 



I 2 
GM ' 



• 2 
JK 



GMr, 



GM 



(3D 



(compare with eqs. 8, 22, and 23). The major change is the equation 
for j'k, which makes the disc dynamics crucially different from the 
relativistic case. 

The specific orbital energy in Newtonian theory is 

•2 



(u r ) 2 

E{r) = ( -^ L + V(r), 



gm 3 - 



(32) 



E — for the parabolic infall above the equatorial plane; the bind- 
ing energy of the disc matter equals —E(r). 

The energy liberated in the disc outside a given radius r equals 
L >r — —E(r)M(r). The total luminosity of the disc is given by 



L = -E(R*)M d , 



(33) 



where M d = M(R*) is the accretion rate through the disc and 
—E(Rt) is the binding energy of matter that reaches the accretor. 
The energy liberated per unit radius is 



Figure 7. The radial velocity of the Newtonian fast disc, u r (r) = dr/dt. 
The long-dashed, dashed, and dotted curves show the cases lo = l cr ~ 
1.204J?*n», lo = i?*M*,and/o = (v / 3/2)i?*n«, respectively (seeeq. 14 
for the definition of l ); u, = (2GM / R*) 1 / 2 . The outer radius of the 
disc, ro = 22q /R*u 2 , is shown by solid dots. The radial velocity of matter 
impinging the disc, u r , is shown by the solid curve (eq. 31). 



dL 
dr 



dr 



[-E(r)M(r)] 



(34) 



The effective potential for radial motion in the disc, V(r), is 
given by equation (32). The important difference from the BH case 
is that the pit at small r is now absent (V — » oo at r — » 0) and the 
potential produces a stronger repulsive centrifugal force. It leads 
to tighter constraints on the angular momentum of the fast disc 
regime. The condition that the disc has no turning point (E > V) 
in terms of angular momentum reads, 



j(r) < j t {r) 



( GM 



+ E 



(35) 



At E = it yields jto = V2jn- The minimal value j t o = R*u* 
at r = i?* implies the necessary condition for the disc formation, 



(36) 



This limitation is twice as stringent as that in the relativistic case 
(compare eq. 36 with eq. 29). 



5.2 Numerical model 

We now compute an illustrative numerical model for the same in- 
flow at Soc as in the relativistic problem, dM/df^oo = M t ot/47r 
and l(8oc) = lo sin#oo. We solve numerically equation (21) cou- 
pled with expressions (17), (20), and (31). The outer radius of the 
disc is given by ro = lo/ GM, and the outer boundary condition 
is u r (ro) = u r (ro). At the inner boundary, r = R*, the accreting 
matter is absorbed by the central object. 

Example solutions are shown in Fig. 7-9. Like the BH case, 
there exists a critical angular momentum, Z cr , such that accretion 
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Figure 8. The streamlines of the Newtonian fast disc, from the outer edge 
onto the surface of the accretor, _R» . The dotted, dashed, and solid curves 
show the cases i = (v / 3/2)il*M* > lo = R*u t , and lo = 'cr ~ 
1.204_R»m*, respectively. The rotation angle, Aip = 2.21 is shown for 
the critical case. 



flow with lo > l CI is stopped by the centrifugal barrier. Numeri- 
cally, we find l CI = 1 .204i?» u* . The minimum angular momentum 
necessary for the disc formation is I, = R*u*/\/2. The fast disc 
regime thus exists at 0.707 < l /R*u* < 1.204. 

At lo — la, the turning point first appears at the inner bound- 
ary of the disc, r = R,. The disc streamlines then touch the accre- 
tor with u r — > 0; the gas moving from the outer edge r rotates by 
an angle A<p « 2.21 (see Fig. 8). The critical behaviour of the BH 
disc is qualitatively different: r cr is twice as large as the BH radius 
and gas rotates an infinite number of revolutions (see Fig. 4). 

The maximum radius of the fast disc is achieved at lo = Z C r, 
r ma x/i?* = l 2 T /GMR„ w 2.90. Hence, the Newtonian fast disc is 
five times smaller than its relativistic counterpart (r max /r 9 « 13.6 
in the BH case). We conclude that the large extension of the BH 
fast disc was caused by purely relativistic effects, namely, by the 
pit in the effective potential. 

We now compute the disc luminosity. First of all, note that the 
total energy released in the process of accretion onto a star of radius 
R, is Ltot = MtotGM / R, . The bulk of the energy is released on 
the star surface, and only a fraction is liberated in the disc itself. 
In Fig. 9 we show the radial distribution of the disc luminosity. In 
contrast to the BH case, the luminosity peaks at r — R* (compare 
with Fig. 5). 

The radiative efficiency of the disc is determined by the bind- 
ing energy, E b = — E(R*) at the inner edge. The E b increases 
with lo and reaches its maximum at lo = l C r, 

1 - 4 ( ^ + " 



E b 



GM 



2/cr 

317 



3/c 



0.0572u;. 




We here accounted for the fact that u r = at r = R* and 
used equation (19). Amusingly, the critical binding energy of the 



Figure 9. The radial distribution of the Newtonian fast disc luminosity 
in units of JWtotu*. The long-dashed, dashed, and dotted curves show the 
cases = 'cr ~ 1.204R»ii*, Iq = -R*m*, and Iq = (v / 3/2)i?*u*, 
respectively. 



Newtonian fast disc in units of «, exactly coincides with the 
binding energy of the marginally stable Keplerian orbit around a 
Schwarzschild BH, (c 2 - E)/c 2 = 1 - 2^2/3 w 0.0572. 



6 DISCUSSION 

6.1 Special features of the fast disc regime 

There are two principle features of the accretion regime studied in 
this paper, which make it different from previous models of rotating 
accretion flows onto black holes: 

(i) The Compton cooled quasi-spherical flow is falling freely 
until it collides with the small-scale disc. The disc is sticky: the 
shock is radiative and pinned to the disc. This regime requires a suf- 
ficiently high accretion rate (see Section 2.3). A similar behaviour 
(held-down shock at high M) was found for spherical accretion 
onto a neutron star (Shapiro & Salpeter 1975). 

(ii) The small-scale disc has a modest angular momentum 
and a high radial velocity. The centrifugal barrier does not prevent 
accretion and matter is swallowed by the black hole in a steady 
regime. The effective throat of the black hole has a large radius of 
2r 3 : the bottleneck for accretion is at 2r g (see Fig. 2). This rela- 
tivistic effect makes the accretion easier and allows the fast disc 
radius to be as large as w 13.6r 9 w 27GM/c 2 (in Newtonian 
gravity it would be five times smaller, see Section 5.2). 

The fast disc is a sticky caustic in a freely falling flow and in 
this respect it resembles the Bondi-Hoyle-Lyttleton (BHL) accre- 
tion (Bondi & Hoyle 1944). In their case, accretion proceeds along 
a one-dimensional caustic called the accretion line. The accretion 
line is fed at each point by matter from an initially plane-parallel 
flow which is axisymmetric and has zero net angular momentum. 
In our case, the accretion disc is fed by an asymptotically isotropic 
inflow with a non-zero net angular momentum. In contrast to the 
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BHL accretion line, the radial momentum of matter impinging the 
disc is directed inwards. As a result, we do not have the degenerate 
stagnation radius that exists in the BHL problem and leads to an 
infinite number of possible solutions. 

The fast disc is crucially different from the standard disc in 
many respects: 

(i) The fast disc forms in accretion flows with a net angular mo- 
mentum l z < 2r g c. 

(ii) The disc is small-scale, ~ 10r 9 . 

(iii) The disc is inviscid and accretes supersonically. 

(iv) The energy release in the disc is caused by the infall absorp- 
tion. 

(v) Most of the energy is released inside 3r 9 . 

(vi) The disc is essentially relativistic; it does not admit any 
Newtonian approximation. 

(vii) The disc has a low mass and its optical depth is not much 
above unity. 

(viii) The disc dynamics depends on the angular momentum dis- 
tribution of the feeding infall. 

6.2 The generated spectrum 

As we argued in Section 1, small-scale discs are likely to form in 
wind-fed X-ray binaries and they are also possible in AGN. We 
now briefly discuss the expected X-ray spectrum emitted by the fast 
small-scale disc. First we note that discs with m > 1 (L > j)Le) 
have optical depths tt > 1. The emission then consists of two 
components: Lh - the hard radiation from the top optically thin 
layer that radiates the vertical kinetic energy of the infall, and L s 
- the radiation produced inside the optically thick disc where the 
horizontal velocities of the infall and the disc get equalised (see 
Section 2.3). Since L s comes out from the optically thick material, 
it is mostly soft radiation. By contrast, Lh is generated by unsat- 
urated Comptonization in the hot layer with the electron tempera- 
ture kT e ~ 100 keV (see eq. 12). The Comptonized X-rays have a 
standard power-law spectrum with a break at ~ 100 keV. The soft 
component L s is the source of seed photons for Comptonization in 
the hot layer. 

The expected two-component spectrum is reminiscent of the 
observed spectra in galactic black hole sources and AGN (see 
Zdziarski 1999 for a review). The observed spectra are well fit- 
ted by two phenomenological models: (i) a hot accretion disc sur- 
rounded by a cold gaseous ring or (ii) a hot corona atop a cold disc. 
Both models are usually considered as modifications of the stan- 
dard a-disc in which the energy is released as a result of viscous 
stresses in a flow with high angular momentum (see Beloborodov 
1999 for a review). The small-scale disc studied here generates a 
similar spectrum and it may provide a preferable explanation of the 
observed spectra in some cases. In general, a Comptonised spec- 
trum is degenerate, i.e. it is insensitive to the heating mechanism 
and/or the accretion dynamics. Detailed studies of the spectrum 
features such as the Fe Kq line and the Compton reflection bump 
from the optically thick disc might help to break the degeneracy. 

6.3 Issues for further study 

(i) Exact calculations of the emitted X-ray spectrum require de- 
tailed modelling of radiative transfer in the post-shock plasma. 

(ii) At high M, the quasi-spherical infall gets so optically thick 
that it will trap the inner radiation and tend to advect it into the BH 
(Begelman 1978). The trapping becomes efficient at m <; 10, when 



the time-scale of radiative diffusion through the infall is longer than 
the accretion time-scale. The corresponding reduction of the ob- 
served luminosity and the impact on the X-ray spectrum should be 
taken into account. 

(iii) Throughout the paper, we assumed that the BH has a small 
spin. The similar problem can be solved for a rotating black hole. 

(iv) We did not consider here the case I > l CT where accretion 
is stopped by the centrifugal barrier. At high /, the fast disc should 
transform into the standard accretion disc. This transition will be 
studied elsewhere (Beloborodov & Illarionov, in preparation). 
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APPENDIX A: CONSERVATION LAWS 

The barion conservation law reads \7k(pu k ) = where V; is the 
covariant derivative in Schwarzschild space-time, p is matter den- 
sity in its rest frame, u l = dx 1 /dX is matter four- velocity, A is 
proper time, and x % are Schwarzschild coordinates. For the flow in 
the thin disc u e = 0, so that u' = {u t ,u r , 0, u v ). The infall has 
4- velocity u l and comoving density p. 

In coordinates x l , the barion conservation law takes the form 
(e.g. Landau & Lifshitz 1975) 

1 d 



(Al) 



y/— g dx k 

where g is the determinant of the metric. The non-zero compo- 
nents of Schwarzschild metric in the equatorial plane are g u = 
-(1 - r g /r)c 2 , g rr = r/(r - r g ), and gee = g vv = r 2 . Then 
sp^g = r 2 . Taking into account that the disc is steady, d/dt = 0, 
and axisymmetric, d/dtp = 0, there are two terms in equation (Al), 
with r— and 6— derivatives. Integrating (Al) over 9 between the 
two surfaces of the disc, one gets 



— — = — 47rr pu 
dr 



M = 2-KrY.u. 



(A2) 



Here S = p^/geH/S.9 is the surface density of the disc (measured 
in its rest frame), A9 = 2H/r is the angular thickness of the disc. 

The energy and momentum conservation reads ViTl — 
where TJ, is the stress-energy tensor. In coordinates x l this equa- 
tion reads (Landau & Lifshitz 1975) 



-g dx k 



I I rpfc\ 1 9gkl rpkl 

( V =-gT i ) = -—T . 



(A3) 



The disc internal energy and pressure are small (heat is radiated 
away) and we neglect their contribution to Tfc. We also neglect vis- 
cous stresses. Tl then takes the simple form 

Tii 2 i 

fe = pc u u k . 

Analogously, the stress-energy tensor of the infall through the disc 
surface is TJ, = pc 2 u t Uk- The stress-energy tensor associated with 
the vertical radiation flux from the disc is Tl = u'qk + Ukq' where 
q* is the four-flux vector (e.g. Misner et al. 1973). Performing ver- 
tical integration across the disc and using equation (A2), we get 



d /,-, \ . dM 2-kF 
(Mm) - Ui— 5-rui 

/ cr 



— — —u u — . 

2 dx l u r 



(A4) 



dr 

Here F — —(dL/dr)(c 2 /ut2nr) is the radiation flux from the 
two faces of the disc (measured in the rest frame of the accreting 
matter). 

For x l = t and x % = ip equation (A4) states the conservation 
of energy, E = —Ut, and angular momentum, j — u v , respec- 
tively, and in both cases the right-hand-side vanishes. The radiative 
term on the left hand-side is ~ L/Mc <§C 1 smaller compared to 
the other two terms and it can be neglected in the conservation law 
of angular momentum. 

For x 1 — r, the right-hand-side in eqiation (A4) is non-zero; 



simple algebra (using utu + u r u r + u v u v 
2r 



-c ) gives 



dgki k i 
——it u — 
uv 



g 2 . (2r-3r 9 ) 2 



r{r-r g ) 



^-u r + 



r 3 (r — r g ) 



Substituting this relation into (A4) and neglecting the radiative 
term, we get equation (21). 



APPENDIX B: THE LUMINOSITY CAPTURE BY THE 
BLACK HOLE 

Assume that each element of the disc dS emits radiation isotrop- 
ically in the rest frame of the accreting matter. Let L c be the lu- 
minosity of dS measured in the rest frame, L\ oc - the luminos- 
ity of dS measured by the local stationary observer (LSO), f3 ~ 
the local disc velocity measured by the LSO (in units of c), and 
7 = (1 — Z? 2 ) -1 / 2 . For a steady disc we have the relation (see e.g. 
Rybicki & Lightman 1979, chapter 4.8) 

d ^ioc _ Lc , R| , 

dQ 4tt 7 4 (1 -p- fi) 3 ' V ' 

where ft is the unit vector specifying direction of an emitted photon 
in the LSO frame and dQ is the element of the solid angle. The local 
luminosity emitted in 47r, as measured by LSO, is equal to the rest 
frame luminosity, L\ oc — L c . 

Consider the local tetrade (e! t ),eLi,ert>,eL)) and let 
03 (r) ,O,/3 (¥,) ) and (ft (r) , Q w , f^>) be the tetrade components 
of (3 and ft, respectively. The /3 (a) (a = 1, 2, 3) is related to the 
four-velocity of the disc, u\ by ft a) = (u Q, /it t )(- 3ciC( /p tt ) 1/2 , 
and 7 = w'(-g tt ) 1/2 . The Q (a) can be written as fi (r) = - cos<9, 
Q w = - sin 6> sin 1,3, n M = sin 9 cos <p, and dfi = Bm0d9d<p. 
Equation (27) now takes the form 



dLi, 



^ioc sin f 



d9d<p 4tt7 4 [1 + cos 9 - @M sin 9 cos <^] 3 

The black hole absorbs radiation emitted within a cone 9 < ^ abs (r) 
(see Misner et al. 1973), where 

27r~ 



§ _ Z7r v (■> r a\ 



6» abs < tt/2 at r > (3/2)r 9 and 6> abs > tt/2 at r < (3/2)r 9 . At 
r = r g , 9 a bs = Ti", i.e. all radiation is absorbed. 

The captured fraction of radiation emitted at a radius r is 



k(Y) = 



2tt 



sin 9 d9 dip 



o 4tt 7 4 [1 + /3M cos 9 - fiM sin 9 cos (pf 



(B2) 



Here /3 (r) = u r /E, (3^ = (cj/rE)(l - r g /r)^ 2 , and 7 = 
(i?/c 2 )(l — r g/ r )^ 1 ^ 2 - The released and captured power from 
dS as measured at infinity is L\ oc (—gtt) and nL\ oc {—gtt), respec- 
tively. The resulting radial distribution of the observed luminosity 
is 



dL obs dL 
-^=d7 [1 - K(r)1 ' 



(B3) 



Equation (B3) is exact for an optically thin disc only. An op- 
tically thick disc partly intercepts its own radiation owing to the 
gravitational bending of light near the black hole, and reemits this 
radiation. As a result of this self-illumination, the radial distribution 
of the disc luminosity changes. We here neglect these effects and 
use equations (B2) and (B3) to estimate the captured luminosity. 
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